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appropriate  results  from  two-dimensional  analyses. 


DST04& 

MELBO  U  R  N  'E 


(C)  COMMONWEALTH  OF  AUSTRALIA  1988 


POSTAL  ADDRESS: 


Director,  Aeronautical  Research  Laboratory, 

P.0.  Box  4331,  Melbourne,  Victoria,  3001,  Australia 


CONTENTS 


Page 

1.  INTRODUCTION .  1 

2.  ANALYTICAL  PROGRAM .  1 

3.  RESULTS .  1 

4.  DISCUSSION  AND  CONCLUSIONS .  2 

5.  REFERENCES .  4 

TABLE 

FIGURES 

DLSTRIBUTION  LIST 
DOCUMENT  CONTROL  DATA 


Acceptor,  ror 

NT!S  CKAI 

01 1<;  i'-n  U 

UlMiViC  .  iCc'<1  G 

Jilitlt.Cd'Ml  _ 

By  .  . 

3  A  ,0'lt 

Av.r.Lbi'ity  r.rvks 


-  1  - 


1.  INTRODUCTION 

In  view  of  the  wide  usage  of  cold-working  and  interference-fits  to  improve  the 
fatigue  performance  of  bolted  joints,  a  programme  is  proceeding  at  these 
laboratories  to  determine  the  stress/strain  fields  which  are  fundamental  to  the 
fatigue  life  enhancement  process. 

Previous  studies  (1,2)  have  examined  the  stress/strain  fields  using  two- 
dimensional  finite-element  analyses  and  also  an  analytical  treatment,  with  particular 
attention  given  to  modelling  interference  by  relative  displacement  loading  (rather 
than  uniform  pressure  loading  which  would  be  correct  only  for  simple  plate 
geometries).  This  paper  extends  the  study  to  a  three-dimensional  finite  element 
analysis  whilst  retaining  relative  displacement  loading. 

2.  ANALYTICAL  PROGRAM 

The  FEM  computations  were  based  on  interference  of  a  high-strength  steel  pin 
in  a  circular  aluminium  plate  ten  times  the  diameter  of  the  pin.  The  plate  thickness 
was  taken  equal  to  hole  diameter.  The  assumed  material  properties  are  as  given  in 
Table  1  and  the  plate  material  has  been  considered  to  have  bilinear  stress/strain 
characteristics,  Fig.  1,  with  isotropic  strain-hardening.  Because  of  the  relatively 
high  yield  point  for  the  pin  material  no  allowance  was  required  for  yielding. 
Interference  levels  of  2%  and  4 %  were  simulated  followed  by  unloading  to  a  cold- 
worked  state  from  both  levels. 

The  FEM  mesh  for  the  plate  was  a  10  degree  sector,  as  shown  in  Fig.  2(a), 
containing  56  twenty-noded,  three-dimensional  brick  elements,  and  18  fifteen-noded 
wedge  elements.  The  pin  mesh,  shown  in  Fig.  2(b),  contained  48  twenty-noded, 
three-dimensional  brick  elements  and  64  fifteen-noded  wedge  elements.  A  small 
hole  (2%  of  the  pin  radius)  was  created  in  the  centre  of  the  pin  mesh  to  avoid  a 
restriction  on  angular  conformation. 

The  first  phase  of  loading,  inducing  interference,  involved  forcing  a  uniform 
relative  radial  displacement  between  pin  and  plate  nodes  at  the  interface  as 
previously  described  in  Ref.  1.  However  representation  of  the  cold-worked  state  was 
most  conveniently  carried  out  by  switching  to  pressure  loading  and  subsequent 
unloading,  using  the  pressure  distribution  in  the  bore  that  was  found  with 
displacement  loading.  Removal  of  the  simulated  interface  pressure  was  equivalent 
to  removing  the  interference. 

The  analysis  used  the  PAFEC  finite  element  scheme  (Level  6)  on  the  ARL 
ELXSI  computer.  The  plasticity  routines  of  that  scheme  are  based  on  Prandtl-Reuss 
equations  in  association  with  the  von  Mises’  yield  criterion. 

3.  RESULTS 

Stress,  strain  and  strain  energy  density  relationships  are  presented  graphically 
along  with  two-dimensional  results  from  Ref.  2  which  are  included  for  comparison. 
The  presentation  of  three-dimensional  results  is  as  follows  : 

Fig.  3.  Circumferential  and  radial  stresses  -  interference. 

Fig.  5.  Circumferential  and  radial  stresses  -  unloaded. 

Fig.  7.  Through-thickness  stresses  -  interference. 
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Fig.  9.  Through-thickness  stresses  -  unloaded. 

Fig.  10.  Stress  distribution  in  bore  -  interference. 

Fig.  11.  Stress  distribution  in  bore  -  unloaded. 

Fig.  12.  Circumferential  and  radial  strains  -  interference. 

Fig.  14.  Circumferential  and  radial  strains  -  unloaded. 

Fig.  16.  Through-thickness  strains  -  interference. 

Fig.  18.  Through-thickness  strains  -  unloaded. 

Fig.  19.  Strain  distribution  in  bore  -  interference. 

Fig.  20.  Strain  distribution  in  bore  -  unloaded. 

Fig.  21.  Strain  energy  at  the  plate  surface. 

Fig.  22.  Strain  energy  at  0.563  of  plate  thickness  from  surface. 

Fig.  23.  Strain  energy  distribution  in  bore. 

Figs.  3,5,7,9,12,14,16  and  18  above  show  the  variation  of  stress  and  strain,  in 
the  radial  direction  at  three  planes  through  the  thickness,  namely,  at  the  surface,  the 
mid-plane,  and  half-way  between.  The  positions  of  mesh  nodes,  including  mid-side 
nodes,  are  shown  on  Fig.  3<a>  to  give  an  indication  of  element  size.  Comparable 
distributions  from  the  two-dimensional  analysis  where  there  is,  of  course,  no 
variation  in  the  thickness  direction,  are  given  as  follows  : 

Fig.  4.  Circumferential  and  radial  stresses  -  interference. 

Fig.  6.  Circumferential  and  radial  stresses  -  unloaded. 

Fig.  8.  Through-thickness  stresses  -  interference  and  unloaded. 

Fig.  13.  Circumferential  and  radial  strains  -  interference. 

Fig.  15.  Circumferential  and  radial  strains  -  unloaded. 

Fig.  17.  Through-thickness  strains  -  interference  and  unloaded. 

It  should  be  observed  that,  in  the  two-dimensional  analyses,  the  mesh  size  in 
the  radial  direction  was  half  that  used  in  the  three-dimensional  analyses  here. 

4.  DISCUSSION  AND  CONCLUSIONS 

(i)  The  circumferential  and  radial  stresses  at  both  2%  and  4%  interference 
in  the  3-D  analysis  show  marked  variations  in  the  thickness  direction 
(Fig.  3).  Although  there  is  little  variation  between  the  values  at  the 
mid-plane  and  the  3/4  thickness  plane,  the  values  at  the  surface  of  the 
plate,  especially  close  to  the  hole,  are  significantly  different.  In 
particular,  at  the  hole,  with  4%  interference  the  radial  stress  at  the 
surface  is  approximately  half  that  at  the  mid-plane.  A  reduction  in  the 
radial  stress  'or  contact  pressure)  at  the  surface  is  to  be  expected, 
because  the  greater  freedom  to  yield  at  the  surface  permits  the 
specified  displacement  to  be  achieved  at  lower  pressure.  The 
occurrence  of  tensile  circumferential  stresses  at  the  surface  idgain 
near  the  hole)  follows  from  the  lower  contact  pressure. 

(ii)  The  stresses  in  the  interior  of  the  plate  (ie.  at  the  mid-  and  3/4- 
thickness  planes)  in  the  3-D  analysis  of  Fig.  3  show  good  agreement 
with  the  values  obtained  in  the  2-D  plane  strain  analysis  of  Fig.  4(b). 

(iii)  The  residual  circumferential  and  radial  stresses  after  load  removal,  as 
given  by  the  3-D  analysis,  are  shown  in  Fig.  5.  Again  there  is  little 
difference  between  the  stresses  in  the  interior  of  the  plate,  but  the 
stresses  at  the  surface  are  noticeably  different;  the  difference  between 
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the  stresses  at  the  surface  and  in  the  interior  is,  however,  less  than  in 
Fig.  3. 

(iv)  The  residual  stresses  in  the  interior  of  the  plate  in  the  3-D  analysis  of 
Fig.  5  are  in  good  agreement  with  the  plane  strain  residual  stresses  of 
Fig.  6(b). 

(v)  In  the  3-D  analysis,  near  the  hole,  neither  the  radial  nor  the 
circumferential  stresses  at  the  surface  change  significantly  between 
2%  and  4%  interference  (Fig.  3).  However,  the  "  waviness”  of  the 
surface  stress  plots  (which  is  even  more  marked  for  the  residual 
stresses  of  Fig.  5)  suggests  that  perhaps  the  finite  element  mesh  size 
near  the  surface  was  not  sufficiently  fine,  especially  in  the  thickness 
direction.  (See  also  (vii)  below.) 

(vi)  The  through-thickness  stresses  (ie.  the  direct  stresses  in  the  thickness 
direction)  are  shown  in  Fig.  7.  Naturally,  these  stresses  are  zero  on  the 
plate  surface.  Under  interference  loading,  the  values  at  the  mid-plane 
are  in  good  agreement  with  the  plane  strain  values  of  Fig.  8(a).  The 
corresponding  residual  stresses,  shown  in  Fig.  9.  exhibit  some 
differences  from  the  plane  strain  values  of  Fig.  8(b);  in  particular,  near 
the  hole  after  2%  cold  working  the  stresses  from  the  3-D  analysis  are 
somewhat  less  than  the  plane  strain  values,  whilst  after  4%  cold 
working  the  3-D  stresses  do  not  exhibit  the  reverse  yielding  apparent  in 
Fig.  8<b). 

(vii)  In  Fig.  10  the  circumferential,  radial,  and  through-thickness  stresses  at 
the  bore  for  interference  loading  are  shown  as  a  function  of  the 
thickness  co-ordinate.  It  can  be  seen  that  all  stresses  exhibit  the  most 
rapid  variation  over  the  element  adjacent  to  the  surface  (the  thickness 
of  which  is  1/16  of  the  plate  thickness.)  It  is  possible  that  use  of  a 
finer  mesh  would  have  modified  both  the  depth  over  which  this  change 
occurred  and  also  the  surface  stress  values.  The  analogous  results  for 
the  residual  stresses  are  shown  in  Fig.  11  and  the  same  trends  are 
evident. 

(viii)  Because  the  finite  element  mesh  size  in  the  radial  direction  is 

relatively  coarse  in  the  3-D  analysis,  the  calculated  locations  of  the 
elastic-plastic  boundaries  may  not  be  accurate.  As  already  remarked, 
in  the  2-D  analysis  a  mesh  size  half  that  of  the  3-D  analysis  was  used. 
Computational  considerations  dictated  the  need  for  a  coarser  mesh  in 
the  3-D  analysis. 

(ix)  Turning  now  to  strains,  it  can  be  seen  from  Fig.  12  that  the 

circumferential  and  radial  strains  at  both  2%  and  4%  interference  in 
the  3-D  analysis  show  much  less  variation  in  the  thickness  direction 
than  do  the  associated  stresses;  in  fact,  the  circumferential  strains 
show  hardly  any  variation  at  all.  (See  also  (xiii)  below.)  This  effect 
may  be  due  in  some  measure  to  the  specified  displacement  applied  at 
the  interface.  The  strains  of  Fig.  12  are  consistent  with  the  plane 
strain  values  of  Fig.  13(b). 
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(x)  Similar  remarks  apply  to  the  residual  circumferential  and  radial 
strains,  those  for  the  3-D  analysis  being  shown  in  Fig.  14,  and  for  the 
plane  strain  analysis  in  Fig.  15(b). 

(xi)  The  through-thickness  strains  (ie.  the  direct  strains  in  the  thickness 
direction  from  the  3-D  analysis  are  shown  in  Fig.  16  for  interference 
loading.  The  very  low  values  found  in  the  interior  of  the  plate  would 
have  been  expected,  being  consistent  with  plane  strain  behaviour. 
However,  the  strains  at  the  surface  are  very  much  less  than  the  2-D 
plane  stress  values  of  Fig.  17(a). 

(xii)  Similar  remarks  apply  to  the  residual  through-thickness  strains,  those 
for  the  3-D  analysis  being  shown  in  Fig.  18,  and  for  the  2-D  plane  stress 
analysis  in  Fig.  17(b). 

(xiii)  In  Fig.  19  the  circumferential,  radial  and  through-thickness  strains  (for 
interference  loading)  at  the  bore  are  shown  as  functions  of  the 
thickness  co-ordinate.  As  already  observed  in  (ix)  above,  the 
circumferential  strain  shows  virtually  no  variation,  although  some 
changes  in  the  other  strains  are  apparent  near  the  surface.  Similar 
remarks  apply  to  the  residual  strains  shown  in  Fig.  20. 

(xiv)  Load  versus  strain  energy  density  plots  (Figs.  21  and  22)  show  a  reverse 
slope  effect  during  the  unloading  phase.  The  possible  significance  of 
this  effect  in  relation  to  fatigue  life  enhancement  has  been  discussed  in 
Ref.  3.  This  effect  disappears  at  a  non-dimensional  radius  of  about  1.5, 
although  that  location  should  be  regarded  as  approximate. 

(xv)  The  strain  energy  density  distribution  at  the  bore  is  shown  as  a  function 
of  the  thickness  co-ordinate  in  Fig.  23.  Markedly  lower  values  are 
apparent  near  the  plate  surface.  This  is  to  be  expected  in  view  of  the 
lower  radial  pressure  near  the  surface,  coupled  with  fairly  constant 
displacements. 
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Table  1.  Assumed  Material  Properties 


Property 

Plate 

Pin 

Modules  of  Elasticity  (MPa) 

69000 

209000 

Poisson’s  Radio 

0.33 

0.30 

Yield  Point  (MPa) 
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1720 

Strain  Hardening  Slope  (MPa> 

1.200 

486 
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FIG.  3.  CIRCUMFERENTIAL  AND  RADIAL  STRESSES  IN  3-D  ANALYSIS  AT 

(a)  2%  INTERFERENCE  AND  (b)  4%  INTERFERENCE.  (MESH  NODE 

POSITIONS  INCLUDING  MID-SIDE  NODES  ARE  SHOWN  ON  FIG.  3(a)) 
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CIRCUMFERENTIAL  AND  RADIAL  STRESSES  IN  2-D  ANALYSIS  A 
2%  AND  4%  INTERFERENCE  ASSUMING  (a)  PLANE  STRESS  AND 
(b)  PLANE  STRAIN  CONDITIONS. 
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FIG.  5.  RESIDUAL  CIRCUMFERENTIAL  AND  RADIAL  STRESSES  IN  3-D 
ANALYSIS  AFTE1'  (a)  2%  COLD-WORKING  AND  (b)  4%  COLD¬ 

WORKING. 
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FIG.  6.  RESIDUAL  CIRCUMFERENTIAL  AND  RADIAL  STRESSES  IN  2-D 
ANALYSIS  AFTER  2%  AND  4%  COLD-WORKING  ASSUMING 
(a)  PLANE  STRESS  AND  (b)  PLANE  STRAIN  CONDITIONS. 
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FIG.  7.  THROUGH-THICKNESS  STRESSES  IN  3-D  ANALYSIS 
(a)  2%  INTERFERENCE  AND  (b)  4%  INTERFERENCE 


s TRESS  STRESS  (KJPb) 


400  o 


FIG.  9 


.o 


-  4  00 


o  A 


-  Surface 

- 3/4  Position 

-  Mid-plane 


-  i  :;oo  o 


'.’o\  dim? nmona:. 


FIG.  9  (a)  AFTER  2%  COLD-WORKING 
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FIG.  9  (b)  AFTER  4 %  COLD-WORKING 

RESIDUAL  THROUGH-THICKNESS  STRESSES  IN  3-D  ANALYSIS 
AFTER  (a)  2%  COLD-WORKING  AND  (b)  4%  COLD-WORKING 
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FIG.  11  (a)  AFTER  2%  COLD-WORKING 


Racial 

i  •  • 

Through- thickness 


-  a  oo 


o 


Circumferential 


-  800 


o 


5  *5 


1 - i - r 


8 


NON  DIMENSIONAL  DEPTH 

FIG.  11  (b)  AFTER  4%  COLD-WORKING 


1  . 


RESIDUAL  STRESSES  ALONG  THE 
(a)  2%  COLD-WORKING  AND  (b) 


BORE  AFTER 
4%  COLD-WORKING 


STRAINS  (n:0-3)  STRAIN  (*10-3) 


D 


Circumferential 


Radial 


—  - Surface 

-  . -  -  3/  4  He  ight 

-  Mid-plane 


FIG.  12  (a)  2%  INTERFERENCE 

Circumferential 


j  Radial 

i  ■  .  ' 

i 

i 

I 

-  Surface 

. —  3/4  Height 

-  Mid-plane 


NON-DIMENSIONAL  RADIUS 
FIG.  12  (b)  4%  INTERFERENCE 


FIG.  12.  CIRCUMFERENTIAL  AND  RADIAL  STRAINS  IN  3-D  ANALYSIS 
AT  (a)  2%  INTERFERENCE  AND  (b)  4%  INTERFERENCE 
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Circumforent ial 


Radial 


NON-DIMENSIONAL  RADIUS 


FIG.  14  (b)  AFTER  4%  COLD-WORKING 


FIG.  14.  RESIDUAL  CIRCUMFERENTIAL  AND  RADIAL  STRAINS  IN  3-D 
ANALYSIS  AFTER  (a)  2%  COLD-WORKING  AND  (b)  4%  COLD 
-WORKING 
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FIG.  15.  RESIDUAL  CIRCUMFERENTIAL  AND  RADIAL  STRAINS  IN  2-D 
ANALYSIS  AFTER  2%  AND  4%  COLD-WORKING  ASSUMING 
(a)  PLANE  STRESS  AND  (b)  PLANE  STRAIN  CONDITIONS 
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FIG.  16  (a)  2%  INTERFERENCE 
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FIG.  16  (b)  4%  INTERFERENCE 


THROUGH-THICKNESS  STRAINS  IN  3-D  ANALYSIS  AT 
(a)  2%  INTERFERENCE  AND  (b)  4%  INTERFERENCE 
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FIG.  18.  RESIDUAL  THROUGH-THICKNESS  STRAINS  IN  3-D  ANALYSIS 
AFTER  (a)  2%  COLD-WORKING  AND  (b)  4%  COLD-WORKING 


Circumferential 
Through- thick ness 

Radial 


NON-  DIMENSIONAL  DEPTH 

19  (a)  2%  INTERFERENCE 

Circumferential 

Through- thickness 


Radial 


a  7  3 

NON -DIMENSIONAL  DEPTH 

19  (b)  4%  INTERFERENCE 


THE  BORE  AT  (a)  2%  INTERFERENCE  AND 


STRAIN  (x  10-3)  STRAIN  (xlO-3) 


'I 

Circumferential 
Through- thickness 


Radial 


9 


non- dimension;*:,  depth 

FIG.  20  (a)  AFTER  2%  COLO-WORKING 
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FIG.  20  (b)  AFTER  4%  COLD-WORKING 

FIG.  20.  RESIDUAL  STRAINS  ALONG  THE  BORE  AFTER 

(a)  2%  COLD-WORKING  AND  (b)  4%  COLD-WORKING 
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FIG.  21.  STRAIN  ENERGY  RESPONSE  AT  VARIOUS  DISTANCES  FROM  THE 
BORE  AT  THE  PLATE  SURFACE.  LOADING  TO  4%  FOLLOWED 
BY  UNLOADING. 
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FIG.  22.  STRAIN  ELL”. GY  RESPONSE  AT  VARIOUS  DISTANCES  FROM  THE 
BORE  ON  A  PLANE  AT  0.563  OF  THE  PLATE  THICKNESS  FROM 
THE  SURFACE  (APPROX  MID-PLANE  POSITION) .  LOADING  TO 
4%  FOLLOWED  BY  UNLOADING. 


DISTRIBUTION 


AUSTRALIA 

Department  of  Defence 

Defence  Central 

Chief  Defence  Scientist 

Assist  Chief  Defence  Scientist,  Operations  (shared  copy) 

Assist  Chief  Defence  Scientist,  Policy  (shared  copy) 

Director,  Departmental  Publications 

Counsellor,  Defence  Science  (London)  (Doc  Data  Sheet  Only) 
Counsellor,  Defence  Science  (Washington)  (Doc  Data  Sheet  Only) 
S.A.  to  Thailand  MRD  (Doc  Data  Sheet  Only) 

S.A.  to  the  DRC  (Kuala  Lumpur)  (Doc  Data  Sheet  Only) 

OIC  TRS,  Defence  Central  Library 
Document  Exchange  Centre,  DISB  (18  copies) 

Joint  Intelligence  Organisation 

Librarian  H  Block,  Victoria  Barracks,  Melbourne 

Director  General  -  Army  Development  (NSO>  (4  copies) 

Aeronautical  Research  Laboratory 
Director 
Library 

Divisional  File  -  Aircraft  Structures 

Author:  R.P.  Carey 

B.C.  Hoskin 

R.  Jones 

G.S.  Jost 

A.  Wong 

M.  Heller 

N. S.  Swansson 

Materials  Research  Laboratory 
Director /Library 

Defence  Science  &  Technology  Organisation  -  Salisbury 
Library 

Navy  Office 

Navy  Scientific  Adviser  (Doc  Data  sheet  only) 

Army  Office 

Scientific  Adviser  -  Army  (Doc  Data  sheet  only) 

Air  Force  Office 

Air  Force  Scientific  Adviser  (Doc  Data  sheet  only) 

Technical  Division  Library 
HQ  Operational  Command  (SMAINTSO) 

HQ  Support  Command  (SLENGO) 


Department  of  Transport  &  Communication 
Library 

Flight  Standards  Division 

Statutory  and  State  Authorities  and  Industry 

Aero-Space  Technologies  Australia,  Manager/Librarian  (2  copies) 

Australian  Airlines,  Library 

Qantas  Airways  Limited 

Ansett  Airlines  of  Australia,  Library 

Hawker  de  Havilland  Aust  Pty  Ltd,  Victoria,  Library 

Hawker  de  Havilland  Aust  Pty  Ltd,  Bankstown,  Library 

Universities  and  Colleges 

Sydney 

G.P.  Steven,  Department  of  Aeronautical  Engineering 
NSW 

Library,  Australian  Defence  Force  Academy 


SPARES  (10  copies) 
TOTAL  (63  copies) 


AL  149 

REVISED  DECEMBER  87 


department  or  defence 


PAGE  CLASSIFICATION 


DOCUMENT  CONTROL  DATA 


UNCLASSIFIED 


PRIVACY  HARKING 


la.  Afl  NUMBER 


[AR-004-579 


It).  ESTABLISHMENT  NUMBER 


ARL-STRUC-TM-478 


2.  DOCUMENT  DATE 


January  1988 


3.  TASK  NUMBER 


pST  87/039 


Three-dimensional  Computation  of 
Stress,  Strain,  and  Strain  Energy 
Density  under  Interference-Fit  and 
after  Cold- Working  of  Holes. 


5.  SECURITY  CLASSIFICATION 

(PLACE  APPROPRIATE  CLASSIFICATION  IN 
BOX  (S)  IE.  SECRET  (SJ.  CONFIDENTIAL  (C) . 
RESTRICTED (RJ .  UNCLASSIFIED  (U) . ] 


6.  NO.  PAGES 

50 


o 


u 


7.  No.  REFS. 

3 


9.  DOWNGRADING/DELIMITING  INSTRUCTIONS 


fR.P.  Carey 


10.  CORPORATE  AUTHOR  AND  AOORESS 


AERONAUTICAL  RESEARCH  LABORATORY 
P.O.  BOX  4331.  MELBOURNE  VIC.  3001 


11.  OFFICE/POSITION  RESPONSIBLE  FOR 

SPONSOR _ “ _ _ 

SECURITY _ 7 _ 

DOWNGRADING  _ 

APPROVAL _ ~ _ 


12.  SECONDARY  DISTRIBUTION  (OF  THIS  DOCUMENT] 


Approved  for  public  release. 


OVERSEAS  ENQUIRIES  OUTSIDE  STATED  LIMITATIONS  SHOULD  BE  REFERRED  THROUGH  ASDIS.  DEFENCE 
DEPARTMENT  OF  DEFENCE.  CAMPBELL  PARK.  CANBERRA.  ACT  2601 

INFORMATION  SERVICES  BRANCH 

13a .  THIS  DOCUMENT  MAY  BE  ANNOUNCED  IN  CATALOGUES  AND  AWARENESS  SERVICES  AVAILABLE  TO... 

No  limitations. 

13b.  CITATION  FOR  OTHER  PURPOSES  (IE.  CASUAL  ANNOUNCEMENT]  MAY  BE  L  1  UNRESTRICTED  OR 

LJ 

j  j  AS  FOR  13a 

14.  DESCRIPTORS 

15.  OROA  SUBJECT 

CATEGORIES 

Finite  element  analysis 

0089D 

Interference  fit  devices 

Fasteners 

Cold  working 

0046E 

Plastic  analysis 

16.  ABSTRACT 


Stress,  strain  and  strain  energy  density  distributions  are  computed  by  three- 
iimensional  finite  element  analysis  for  interference  fit  and  cold-worked 
conditions  of  a  steel  pin  in  an  annular  aluminium  alloy  pic  te,  the  plate  thickness 
oeing  equal  to  pin  diameter.  Comparisons  are  made  with  the  appropriate 
results  from  two-dimensional  analyses. 


THIS  PAGE  IS  TO  BE  USED  TO  RECORD  INFORMATION  WHICH  IS  REGl/IREO  BY  THE  ESTABLISHMENT  FOR 
ITS  OWN  USE  BUT  WHICH  WILL  NOT  BE  ADDED  TO  THE  OISTIS  DATA  UNLESS  SPECIFICALLY  REQUESTED. 


16.  ABSTRACT  (COffT.) 


AERONAUTICAL  RESEARCH  LABORATORY,  MELBOURNE 


18.  OOCUH6KT  SERIES  ANO  NUMBER 

Aircraft  Structures 
Technical  Memorandum  478 


19.  COST  cooe 

24  1160 


20.  TYPE  OF  REPORT  ANO  PER I 00  COVERED 


21.  COMPUTER  PROGRAMS  USED 

PAFEC 


22.  ESTABLISHMENT  FILE  REF.  (SI 


23.  additional  Information  (as  required) 


